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ABSTRACT 



o 
u 

O . We review the results of high dispersion spectroscopy of 17 circumnuclear star- 

-(_> ' forming regions (CNSFRs) in 3 nearby early spiral galaxies, NGC 2903, NGC 3310 and 

^ ; NGC 3351. 

Wc find that single Gaussian fitting to the H/3 and [Oiii]A5007A line profiles 
results in velocity dispersions around 32knis~^ and 52knis~^, respectively, while the 
IR Can triplet cross-correlation technique provides stellar velocity dispersion values 
^ • close to 50kms~^. 

^^ \ Even though multiple kinematical components are present, the relation between 

gas velocity dispersion and Balmcr emission line luminosity (L-cr relation) reproduces 
the correlation for disk giant Hii regions albeit with a larger scatter. The scatter in 
the L-(T relation is considerably reduced when theoretical evolutionary corrections are 
'^ ■ applied suggesting that an age range is present in the sample of CNSFRs. 

^^ ' To analyse the observed complex profiles, we performed multiple Gaussian compo- 

nent fits to the H/3 and [Oiii]A5007A lines obtaining optimal fits with two Gaussians 
of different width. These best fits indicate that the narrower component has aver- 
age velocity dispersion close to 23kms~"'^ while the broader component shows average 
values in the range 50-60 km s""'^ for both lines, close to the observed stellar velocity 
dispersions. The fiuxes of the broad and narrow H/3 components are similar. This is 
not the case for [Oiii]A5007 A for which the broad components have higher fluxes than 
the narrow ones, thus producing a clear segregation in their [Olll]/H/3 ratios. 

We suggest a possible scenario for understanding the behaviour of CNSFRs in the 
L-cr and cFgas-c* diagrams involving an inner gaseous disk responsible for the narrow 
component of the emission lines. 

Our main conclusion is that the presence of different kinematical components with 
similar total fluxes in the emission line spectrum of CNSFRs raises important doubts 
regarding properties of the ionized gas derived from global line ratios obtained with 
low resolution spectroscopy in star-forming regions in the central regions of early type 
galaxies. Given the ubiquity of these star-forming systems, ionized gas analyses should 
be done preferably from high dispersion spectra with high spatial resolution. 

Key words: HII regions - galaxies: kinematics and dynamics - galaxies: starburst - 
galaxies: star clusters. 



1 INTRODUCTION 
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active nucleus. These structures are frequently arranged in a 
ring pattern with a diameter of about 1 kpc (Morgan 1958; 
Sersic & Pastoriza 1965, 1967) and their associated star for- 
mation rate (SFR) is considerably higher than the average 
observed in galactic discs. In fact, in the ultraviolet (UV), 
the luminosity coming from their young massive stars can 
dominate the observed circumnuclear emission even in the 
presence of an active nucleus (Gonzalez-Delgado et al. 1998; 
Colina et al. 2002). 

The properties of circumnuclear star-forming regions 
(CNSFRs) are different enough from giant extragalactic Hll 
regions that they can be considered to constitute a distinc- 
tive mode of star formation in galaxies (Kennicutt 1998). 
This distinct nature was fully revealed with the opening of 
the mid and far-IR spectral ranges (see for example Rieke 
& Low 1972; Harper & Low 1973; Rieke & Lebofsky 1978; 
Telesco & Harper 1980). CNSFRs look more compact and 
show higher peak surface brightness than disc Hll regions 
(Kennicutt et al. 1989) while their optical spectra show the 
weak oxygen forbidden lines characteristic of a high abun- 
dance (Diaz et al. 2007, hereafter D07). 

It has been suggested that the substructures (sizes of 
the order of 10 pc or less) found in circumnuclear star form- 
ing complexes correspond to massive (M ~ 10^ Mq) young 
(age < 1 Gyr) clusters. If this is the case, they can be mod- 
elled as single stellar population applying population syn- 
thesis techniques and their properties such as age, chemical 
evolution and star formation histories could be derived with 
only small uncertainties due to discreetness related scatter 
(see e.g. Cerviiio et al. 2002) . Although there have been 
some recent works devoted to the study of CNSFRs most 
of them have concentrated in deriving the SFR and ages 
of these structures and little has been done with respect to 
their structure and internal kinematics. 

In the last few years we have obtained spectroscopic 
data of CNSFRs in nearby galaxies with a resolution enough 
to allow the measurement of stellar and gaseous velocity dis- 
persions that, combined with their sizes obtained from the 
Hubble Space Telescope (HST) imaging, have yielded val- 
ues for the dynamical masses of these objects (Hagele et al. 
2007, 2009, 2010, hereafter H07, H09 and HIO, respectively). 
In this paper we present the results obtained from the joint 
analysis of the data and address several aspects related to 
the properties of these structures. In particular, we analyse 
the relations between stellar and gaseous velocity dispersions 
and the behaviour of the velocity dispersion-luminosity (cr- 
L) and mass-luminosity (M-L) relations for CNSFRs. We 
also draw attention to several concerns related to conclu- 
sions about gas properties of central regions of galaxies, like 
average chemical abundances or characteristics of ionizing 
sources, derived from data with low spatial and spectral res- 
olution. 



2 SUMMARY OF DATA AND DERIVED 
PARAMETERS 

The CNSFRs studied in the present work correspond to 
three early type barred spiral galaxies, NGC 2903, NGC 3310 
and NGC 3351 whose properties are given in Table 1. 

For 17 star forming complexes we obtained high res- 
olution spectroscopy in the blue (4779 to 5199 A) and red 



Table 1. The galaxy sample. 



Property 


NGC 2903 


NGC 3310 


NGC 3351 


RA (2000)'' 


09 32 10.1 


10 38 45.9 


10 43 57.7 


Dec. (2000)" 


-1-21 30 03 


+53 30 12 


+11 42 14 


Morph. type 


SBbc 


SABbc 


SBb 


Distance (Mpc) 


8.6 ±0.5'' 


15 + 1" 


10.05 ±0.88^= 


pc/arcscc 


42 ±2 


73 ±5 


49 ±4 


Bt (mag)° 


9.7±0.1 


11.2±0.1 


10.5 ±0.1 


E(B-V)g,,(mag)'' 


0.031 ±0.01 


0.030 ±0.01 


0.028 ±0.01 



" dc Vaucoulcurs ct al. (1991) 
*> Bottinelli ct al. (1984) 
•^ Graham et al. (1997) 



(8363 to 8763 A) with resolutions of 0.21 and 0.39 A per 
pixel, respectively, providing a comparable velocity resolu- 
tion of about 13kms~'^. The blue spectral region comprises 
the emission lines of H/3A486lA and [Olll]A5007A and 
have been used to measure the gas velocity dispersion. The 
red spectral region comprises the stellar absorption lines of 
the Can triplet (CaT) at AA8494, 8542, 8662 A, and have 
been used to measure the stellar velocity dispersion. The two 
spectra have been obtained simultaneously using the double 
arm Intermediate-dispersion Spectrograph and Imaging Sys- 
tem (ISIS) attached to the William Herschel Telescope and 
therefore we are confident that they refer to the same re- 
gions. More details about the observations can be found in 
our previous works (H07, H09, HIO). 



The velocity dispersions of the gas, Og 
lated as 



were calcu- 



^ - ^f (1) 

where Gm and ai are the measured and instrumental disper- 
sions, respectively, in kms~^. Oi was measured directly from 
the sky emission lines and is about 10.4 km s~^ at A 4861 A. 
No thermal broadening correction was applied, thus our sig- 
mas slightly overestimate the gas intrinsic dispersion. 

The single component Gaussian fits performed to the 
H/3 line were not totally satisfactory. In fact, in all cases, 
the optimal fit was found for two different Gaussian com- 
ponents. In most cases, similar components also provided 
optimal fits to the [Olll] lines. Fig. 1 shows fits to the H/? 
lines in CNSFRs of each of the three galaxies studied as 
illustrative examples. 

The stellar velocity dispersions, a», were estimated from 
the CaT lines by cross-correlating the spectrum correspond- 
ing to each of the regions with stellar templates. An example 
of the CaT spectral region can be seen in Fig. 2. 

When observed with the HST resolution, most of the 
observed CNSFRs seem to be a composite of more than one 
cluster with radii between 1.5 and 6.2 pc. The virial theorem 
has been applied to estimate upper limits to the dynamical 
masses inside the half light radius for each identified knot 
under the assumption that the clusters are gravitationally 
bound, spherically symmetric and have isotropic velocity 
distributions [cr^ (total) ^ 3 al\. The general expres- 
sion for the virial mass of a cluster is 



Kh = v^IR/G 



(2) 



where R is the effective gravitational radius and rj is 
a dimensionless number that takes into account de- 
partures from isotropy in the velocity distribution 
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Figure 1. Optimal Gaussian fits to tlie H/3 line for observed regions in the galaxies: NGC 2903, NGC 3351 y NGC 3310. The two 
different components labelled "broad" and "narrow" are sfiown. 
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Figure 2. Spectrum of region R4 in NGC 2903 in the IR Call 
triplet spectral region. 



and the spatial mass distribution, binary fraction, 
mean surface density, etc. (Boily et al. 2005; Fleck 
et al. 2006). Following Ho &: Filippenko (1996a,b), 
and for consistence with HOT, H09, HIO and Hagele 
(2008), we obtain the dynamical masses inside the 
half-light radius using 77 = 3 and adopting the half- 
light radius as a reasonable approximation of the 
effective radius. Since our measured stellar velocity dis- 
persion encompass each whole region to which the single 
clusters belong, their masses can be overestimated to some 
extent. However we think the mass overestimate cannot be 
large, because the equivalent widths of the CaT lines point 
to the spectral energy distribution of each region being dom- 
inated by the young stellar population belonging to the clus- 
ters rather than by the galaxy bulge. The adopted mass for 
each region has been obtained by summing up the individ- 
ual masses of the single clusters. More details about these 
procedures can be found in H07. 

Table 2 lists the velocity dispersions in each CNSFR 
together with the derived values of the dynamical masses. 
Column 1 gives the name of the region as coded in the orig- 
inal publication; colunrn 2 gives the stellar velocity disper- 
sion {(T*); column 3 gives the value of the gas velocity dis- 
persion derived using a single Gaussian fit [a a) to the H/3 
line; columns 4 and 5 give the corresponding values derived 
by the two-component optimal fit that have been labelled 
narrow ((t„) and broad ((Jb), respectively. Columns 6 to 8 



provide the same information for the [Olll] emission lines. 
Finally, column 9 gives the derived dynamical mass. Mast- 

There are three objects with two sets of values listed 
in the table which correspond to different observations ac- 
quired in two different slit position angles. 

For every region we have derived the mass of the 
ionising stellar population, (Afion), from published values 
of the Ha luminosity (Planesas et al. 1997; Diaz et al. 
2000; Pastoriza et al. 1993) corrected for reddening as- 
suming the galactic extinction law of Miller & Mathews 
(1972) with Rv = 3.2 and, when necessary, correcting for 
recent changes in their distances. The solar metallicity 
[12-f log(0/H) = 8.919] single burst models by Garcia- Vargas 
et al. (1995) were used. These models assume a Salpeter 
initial mass function (Salpeter 1955, IMF) with lower and 
upper mass limits of 0.8 and 120 Mq and provide the num- 
ber of ionising Lyman continuum photons [Q(^o)] per unit 
mass of the ionising population, [Q{Ho)/Mion\. This number 
decreases with the time evolution of the population, and is 
related to the computed equivalent width of the H/3 line as 
(e.g. Diaz et al. 2000): 



log [QiHo)/M^on] = 44.48 + 0.86 log [EW{RI3)] 



(3) 



The total number of ionising photons for a given re- 
gion has been derived from the Ha luminosities (Leitherer 
& Heckman 1995): 



Q{Ho) = 7.35 



10" L(Ha) 



(4) 



The final expression for the derivation of Mion is then: 



M,on 



7.35 X 10"L(Ha) 

]^Q44.48 + 0.86 log[BW(H/3)] 



(5) 



The equivalent widths of the H/3 line have been taken 
from Pastoriza et al. (1993) and D07. Given that the H/3 
equivalent widths are probably affected by an underlying 
older stellar continuum not belonging to the ionizing cluster 
itself, the listed masses for these clusters should be consid- 
ered upper limits. On the other hand, no photon escape or 
absorption of ionizing photons by dust has been considered, 
which would place lower limits to the mass. These two effects 
probably cancel each other to a certain degree. 

The mass of ionised gas (Mhii) associated to each star- 
forming region complex was derived from the Ha luminosi- 



4 G. F. Hdgele et al. 



Table 2. Velocity dispersions and dynamieal masses for the observed CNSFR 









H/3 






Oiii] 






Region 


0-, 


O-s 


CT,i 


o-fc 


CTs 


CTn 


o-fc 


AL 










NGC2903 








R1+R2 


60±3 


34±2 


23±2 


51±3 


73 ± 8 


26 ± 8 


93 ± 9 


1816±48 


R1+R2 


64±3 


35±2 


27±2 


53±4 


71 ± 9 


27 ± 7 


95 ± 10 


2054±45 


R4 


44±3 


32±2 


20±2 


47±4 


76 ± 10 


35 ± 9 


89 ± 8 


853±28 


R7 


37±3 


32±4 


29±5 


34±8 


59 ± 10 


17 ± 5 


67 ± 8 


642±26 










NGC3310 








R1+R2 


80: 


33±4 


24±5 


54±7 


31 ±4 


22 ± 5 


50 ± 4 


911: 


R4 


36±3 


34±3 


28±4 


55±5 


32 ± 3 


26 ± 3 


52 ± 4 


886 ±30 


R4+R5 


38±3 


27±4 


22±4 


46±5 


22 ± 3 


18 ± 3 


40 ± 2 


1317 ±41 


R6 


35±5 


30±3 


23±3 


54±7 


28 ± 3 


21 ± 3 


56 ± 3 


397 ±33 


S6 


31±4 


27±3 


20±3 


47±4 


26 ± 3 


19 ± 3 


47 ± 3 


210 ±17 


R7 


44±5 


21±5 


18±4 


41±4 


17 ±4 


14 ± 3 


36 ± 3 


1413±60 


RIO 


39±3 


38±3 


26±2 


54±2 


40 ± 3 


26 ± 3 


59 ± 4 


588 ±28 


J 


— 


34±2 


25±2 


61±3 


30 ± 2 


22 ± 2 


57 ± 3 


— 










NGC3351 








R2 


50±1 


26±1 


17±3 


45±3 


72 ± 7 


21 ± 4 


74 ± 5 


129±6 


R2 


51±6 


29±3 


16±2 


43±2 


69 ± 9 


23 ± 5 


76 ± 8 


131±23 


R3 


55±5 


35±1 


25±3 


59±4 


67 ± 7 


28 ± 4 


71 ±4 


417±31 


R3 


59±7 


39±5 


24±3 


59±3 


70 ± 7 


24 ± 6 


74 ± 9 


477±47 


R4 


66±4 


37±4 


29±3 


65±4 


76 ± 8 


— 


— 


87±12 


R5 


47±4 


34±2 


30: 


76: 


56 ± 7 


— 


— 


49±8 


R6 


39±6 


29±6 


16±3 


46±4 


46 ± 7 


— 


— 


434±39 



Note: velocity dispersions in kms~ 
corresponding errors indicate that 
errors (greater than 40%). 



-^, masses in 10^ Mq . The colons instead of the 
those quantities have large associated 



ties, using the electron density (Ne) dependency relation 
given by Macchetto et al. (1990) for an electron tempera- 
ture of 10* K. 



Afmi = 3.32 X 10"^" L(Ha) iV; 



(6) 



The electron density for each region was taken from Pas- 
toriza et al. (1993) and D07. For those regions of NGC3310 
without electron density estimates a value of N^ equal to 
100 cm~ , corresponding to the average value of other CNS- 
FRs of this galaxy, was assumed. 

The values of the Hq luminosities for each CNSFR to- 
gether with quantities derived from them are given in Table 
3. Also in this table, and following the work by Relafio & 
Beckman (2005), we list the individual Gaussian component 
fluxes as fractional emission (EM/), in percentage, relative 
to the total line flux. 

These have been calculated from the emission line fluxes 
estimated for each H/3 component from the parameters de- 
rived using the fitting from the NGAUSS task of IRAF^. Indi- 
vidual component Ha luminosities follow directly from those 
with the only assumptions that the reddening corrections are 
the same for both components (although this is not neces- 
sarily true, see Hagele et al. 2012) and that the EM/ remains 
constant across the whole nebula. It must be noticed that 
the slit is 1 arcsec wide and the estimated diameter of the 
entire star-forming complexes from which the photometric 
fluxes were estimated, are about 2 arcsec. 

As in the case of Table 2 there are several objects with 



^ Image Reduction and Analysis Facility, distributed by NOAO, 
operated by AURA, Inc., under agreement with NSF. 



two sets of values listed in Table 3 which correspond to two 
different observations as explained above. 



3 DISCUSSION 

3.1 Star and gas kinematics 

Given the emission line profile differences between permitted 
and forbidden lines already reported by us in previous works 
(H07, H09, HIO), it is important to study the relationship 
between the velocity dispersion of the stars and the ionised 
gas for all the CNSFRs that have been observed at high 
dispersion and with similar instrumentation. 

A quick inspection of Table 2 shows that while the stel- 
lar velocity dispersion in the CNSFRs ranges between ~30 
and ~75 kms~^ with an average value of 49 kms~^, the ion- 
ized gas velocity dispersion, as obtained from the H/3 emis- 
sion line using a single Gaussian fit, shows a range between 
~20 and ~40 kms"'^ with an average value of 32 kms"'^. On 
the other hand, the [Olll] velocity dispersion, also measured 
from a single Gaussian fit, that ranges between ~17 and 
~75 kms"'^ with an average value of 52 kms~^ has values 
closer to the stellar velocity dispersion perhaps indicating 
some extra source of broadening with respect to H/? (see a 
detailed discussion in H09). 

The systematic shift between the ionized gas and stellar 
velocity dispersion is clearly seen in Fig. 3 that shows the 
gas velocity dispersion measured by fitting a single Gaus- 
sian to the H/3 and [Olll] lines versus the stellar velocity 
dispersion for all the regions. Even though the scatter is 
large it is possible to see that while most of the CNSFRs 
with stellar velocity dispersion below --^40 kms^^ (all but 
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Table 3. Ha luminosities and derived quantities, and H/3 equivalent widths. 



Region 


L{na)t 


EM/„ 


EM/i, 


L(na)n 


L{na)t 


Mio„ 


Mhii 


EW(H^) 










NGC 2903 










R1+R2 


66.3'' 


49 


51 


32.3 


34.0 


18.9 


0.79 


12.1 


R1+R2 


— 


57 


43 


38.1 


28.2 


— 


— 


12.1 


R4 


38.9" 


32 


68 


12.6 


26.3 


24.6 


0.48 


4.8 


R7 


31.3" 


59 


41 


18.6 
NGC 3310 


12.7 


23.6 


0.30 


3.9 


R1+R2 


102'' 


48 


52 


49.3 


52.7 


13.9 


3.38 


28.6 


R4 


144* 


58 


42 


83.5 


60.5 


17.6 


4.78 


32.4 


R4+R5 


218'' 


62 


38 


136 


82.0 


21.4 


7.24 


41.7 


R6 


57.3'' 


53 


47 


30.5 


26.8 


12.4 


1.90 


16.7 


S6 


62. S'^ 


54 


46 


34.0 


28.5 


17.4 


2.07 


12.5 


R7 


45.5'' 


72 


28 


32.7 


12.8 


8.7 


1.51 


19.4 


RIO 


45.5'' 


39 


61 


18.0 


27.5 


15.7 


1.51 


9.7 


J 


573'' 


51 


49 


294 
NGC 3351 


279 


31.4 


9.52 


82.5 


R2 


28.3" 


40 


60 


11.4 


16.9 


9.93 


0.21 


9.5 


R2 


— 


35 


65 


10.0 


18.3 


— 


— 


9.5 


R3 


70.0" 


49 


51 


34.4 


35.7 


15.3 


0.54 


16.5 


R3 


— 


33 


67 


22.9 


47.1 


— 


— 


16.5 


R4 


23.2" 


48 


52 


11.1 


12.1 


6.23 


0.25 


13.0 


R5 


10.3" 


80 


20 


8.2 


2.1 


6.17 


0.09 


5.1 


R6 


7.5" 


46 


54 


3.4 


4.0 


8.88 


0.07 


2.3 



Note: EM/ in percentage, luminosities in 10'^** ergs"'^, masses in 10^ Mq. 

"From Planesas ot al. (1997) corrected for the different adopted distances, and for 

reddening using E(B-V) from Perez-Olea (1996). 

*> From Diaz et al. (2000). 

" From Pastoriza et al. (1993). 



two belonging to NGC 3310) show similar stellar, H/3 and 
[Olll] velocity dispersions, the CNSFRs with stellar velocity 
dispersion above ~40 kms^^ show a clear dichotomy with 
the [Qui] profiles being systematically broader than those 
of H/3. The Balmer emission lines are affected by the 
absorption produced by the underlying stellar pop- 
ulation that depress the lines (see a detailed dis- 
cussion about this effect in H07, H09 and HIO). To 
minimise the errors introduced by the underlying 
population we defined a pseudo-continuum at the 
base of the line to fit the emission lines (for details, 
see Hagele et al. 2006). Nevertheless, the presence 
of a conspicuous underlying population could still 
be affecting the agas derived from the H/3 emission 
lines. On the other hand, the fact that the regions 
with higher values of cr, have systematically larger 
disagreement with the agas ^= a, line could be related 
with that they are expected to be the most massive 
ones having the larger amount of stars and the most 
conspicuous underlying stellar populations. 

As mentioned above, most of the observed CNSFR 
emission line profiles show wings detectable in both the H/? 
hydrogen recombination line and the [Olll] forbidden line. 
These wings can be attributed to the presence of a compo- 
nent broader than the core of the line. To investigate the 
origin of this broad component we have assumed that it can 
be represented by a Gaussian profile. It should be kept in 
mind, however, that given that we do not have an a-priori 
scenario for the origin of the broad component, there is no 
justification to assume a Gaussian profile either, apart from 
the fact that this assumption simplifies the analysis. If we 



assume that the wings are due to stellar winds, the broad 
component profile should be that of an expanding shell, i.e. 
a very flat top profile in which case the fluxes of the broad 
component obtained with a broad Gaussian will be overes- 
timated and the dispersion of the narrow component will 
be underestimated. Yet, the best Gaussian fits involved two 
different components for the gas: a narrow component with 
velocity dispersions lower than the stellar one and a broad 
component with a velocity dispersion similar to or slightly 
larger than that measured for the stars. The resulting values 
of the multiple Gaussian fit to the observed profiles for both 
H/3 and [Olll] narrow and broad components are listed in 
Table 2 and Table 3. 

In Fig. 4 and Fig. 5 we show the velocity dispersion 
of the narrow and broad Gaussian components of H/3 and 
[Qui], respectively, versus the stellar velocity dispersion. 
The narrow components of both H/3 and [Qui] span a sim- 
ilar range in velocity dispersion with an average value of 
~23 kms^^. On the other hand, the values of the velocity 
dispersion of the broad component of H/3 are similar to those 
of the stellar velocity dispersion while the broad component 
of the [Olll] Une shows, in most cases, values of a in excess 
of the stellar one. The correlation coefficient of the broad 
component of both H/3 and [Olll] with the stellar velocity 
dispersion is much smaller than 0.5 indicating basically no 
correlation. 

Bearing in mind that we may over-interpret the results 
of what is in reality a small sample, the fact that the narrow 
component of the gas in the CNSFRs in 3 different galaxies 
is, inside the error bars, roughly constant suggests that it 
may be originated in an inner disc that would be mainly 
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Stellar o., vs Gas 0(H|i) , (^([oiiii] total 




Stellar o, vs o, 



narrow and broad 




Log (a.) 



Log (a.) 



Figure 3. Gaseous versus stellar velocity dispersions from single 
Gaussian fits. Pluses correspond to H/3 while open diamonds cor- 
respond to [Qui]. The continuous line represents: cTgas = cf stars ■ 



Figure 5. Gaseous (0[III]) versus stellar velocity dispersions for 
all the observed circumnuclear regions. The continuous line rep- 



resents: (T„ 




Log (tj.) 

Figure 4. Gaseous (H/3) versus stellar velocity dispersions for all 
the observed circumnuclear regions. The continuous line repre- 
sents: (Tnas -— ^stars' 



supported by rotation. On the other hand, the stars and the 
gas responsible for the broad component could be related to 
the star-forming regions themselves and therefore would be 
mostly supported by dynamical pressure (see Pizzella et al. 
2004, and references therein). This could, in principle, cor- 
respond to the gas response to the gravitational potential 
of the stellar cluster, thus explaining the coincidence with 
the stellar velocity dispersion in most cases for the broad 
component of the H/3 line. To disentangle the origin of the 



two different kinematical components it will be necessary to 
map these regions with higher spectral and spatial resolu- 
tion and much better signal-to-noise ratio in particular for 
the [Olll] lines. 

The development of the Integral Field Unit 
(IFU) instruments to perform 3D spectroscopy dur- 
ing the last years has provided the spatial cov- 
erage required to study extended galactic or ex- 
tragalactic star-forming regions (see e.g. Monreal- 
Ibero et al. 2010, 2011; Rosales-Ortega et al. 2010, 
2011; Garci'a-Benito et al. 2010; Sanchez et al. 2011; 
Lopez-Sanchez et al. 2011; Perez-Montero et al. 
2011; Kehrig et al. 2012; Sanchez et al. 2012; Lopez- 
Hernandez et al. 2013; James et al. 2013; Firpo et al. 
2013). Recently the first results of very powerful sur- 
veys of galaxies using IFU instruments have been 
released. The Calar Alto Legacy Integral Field Area 
(CALIFA) survey is obtaining spatially resolved 
spectroscopic information of a selected sample of 
~600 galaxies in the Local Universe (Sanchez et al. 
2012). This survey has an effective spatial resolution 
of about 2 arcsec per fibre, and the final spectral res- 
olution in FWHM is ~6.5 A, (i.e. a ~ 150 kms"^) for 
the V500 grating and '^2.7A, (i.e. a ~ 85 kms"^) 
for the V1200. These instrumental configurations 
are not enough to disentangle the origin of the two 
different kinematical components w^ith the required 
accuracy. The possibility of a future survey per- 
formed using the Sydney-AAO Multi-object Inte- 
gral field spectrograph (SAMI) w^ould be better for 
our purposes since its higher spectral resolution is 
R=A/AA~ 13000 (i.e. a ~ 23 kms"^ at 4861 A), and 
each fibre subtends 1.6 arcsec on the sky (Croom 
et al. 2012). How^ever, medium or high dispersion 
slit spectroscopy are a better option for spectropho- 
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Table 4. Logarithm of the [Olll]5007/H/3 line ratio. 



NaiTow component 



Region 


log([0 


n]5007/H/3) 






One eomponent 


Narrow 


Broad 




NGC 2903 




R1+R2 


-1.03±0.05 


-1.64±0.09 


-0.76±0.10 


R1+R2 


-0.92±0.07 


-1.51±0.11 


-0.57±0.14 


R4 


-0.78±0.11 


-l.OlitO.lO 


-0.71±0.15 


R7 


-0.78±0.10 


-1.79±0.12 


-0.40±0.18 




NGC 3310 




R1+R2 


0.42±0.01 


0.39±0.01 


0.45±0.03 


R4 


0.28±0.01 


0.28±0.01 


0.28±0.04 


R4+R5 


0.42±0.01 


0.39±0.01 


0.44±0.03 


R6 


0.28±0.01 


0.30±0.01 


0.28±0.06 


S6 


0.21±0.01 


0.21±0.01 


0.22±0.08 


R7 


0.23±0.01 


0.21±0.01 


0.32±0.09 


RIO 


0.19±0.01 


0.20±0.04 


0.20±0.08 




NGC 3351 




R2 


-1.07±0.06 


-1.66±0.08 


-0.93±0.12 


R2 


-1.01±0.06 


-1.55±0.08 


-0.96±0.13 


R3 


-I.10it0.06 


-1.57±0.07 


-0.93±0.10 


R3 


-1.00±0.06 


-1.52±0.09 


-0.89±0.10 



tometry, or when the object is very compact, or even 
extended but with few star-forming knots. This is 
also the case w^hen good spatial and spectral reso- 
lution and simultaneous wide spectral coverage are 
required (see e.g. Gumming et al. 2008; Firpo et al. 
2010, 2011; Hagele et al. 2006, 2007, 2008; Hagele 
2008; Hagele et al. 2009, 2010, 2011, 2012; Lopez- 
Sanchez & Esteban 2009, 2010a, b; Lopez-Sanchez 
2010; Amorin et al. 2012). 
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Figure 6. Histogram of the logarithmic [Oin]/H/3 ratio for the 
different gaseous kinematie eomponents compared to values de- 
rived from single Gaussian fits. 



3.2 Effects on abundance estimates 

The presence of two kinematically segregated gaseous com- 
ponents could have an effect on the classification of the 
activity in the central regions of galaxies through diag- 
nostic diagrams and also on abundance determinations if 
they are made from low dispersion spectra by biasing the 
line ratios involved (see e.g. Esteban & Vilchez 1992; 
Lopez-Sanchez et al. 2007; Hagele et al. 2012; James 
et al. 2013,?). The only line ratio available from our 
data for the two kinematically separated components is 
[Olll]/H/3, whose logarithmic ratios are listed in Table 
4. We find that in all the CNSFRs of NGC 2903 and 
NGC 3351, ([Oiii]/H/?)„a„.<,^. < ([Oiii]/H/3)6,o,d which indi- 
cates a higher excitation for the broad component gas. In the 
case of NGC 3351 the mean value of log ([Olll]/H/3)„arroii) 
is -1.57 as compared to -0.93 for the same line ratio of the 
broad components, while a value of -1.04 is obtained from a 
single Gaussian fit (see Fig. 12 in H07). In these cases the 
broad component contributes more than 80% to the [Olll] 
total flux while this contribution is smaller than 40% for the 
H/3 Une. This is not the case for the CNSFRs of NGC 3310 
for which the contribution by the broad component to both 
the [Olll] and H/3 emission line fluxes is similar and about 
50% which maintains the [Olll]/H/3 ratio almost constant. 

Figure 6 shows a histogram of the logarithmic [Olll]/H/? 
ratio for the different components of the gas, broad and nar- 
row, and also for the case of a single Gaussian fit (one com- 



ponent). It can be seen that for the CNSFRs of NGC 2903 
and NGC 3351, the gas showing different velocity dispersions 
also show different ionization properties with the broad com- 
ponent one showing a higher degree of excitation (see also 
Table 4). The data with the highest [Olll]/H/3 ratio corre- 
spond to the regions in NGC 3310 for which no differences 
are found among the different cases. 



3.3 cr-L and M-L relations 

Fig. 7 shows the distribution of the CNSFRs in the L(Ha) 
vs. (7 plane. In the left column, upper panel we plot the log- 
arithm of the total Ha luminosity (values taken from the 
literature corrected for extinction) versus the gas velocity 
dispersion derived using a single Gaussian fit. In the mid- 
dle and lower panels we show the same relation but for the 
Hq luminosity of the narrow [Ln(HQ)] and broad [L6(Hq)] 
components of the two-Gaussian fits derived using the EM/ 
estimated from the fit to the H/3 emission line, versus the cor- 
responding velocity dispersion given by the multi- Gaussian 
fit (an and ab, respectively). The linear fits from Bosch et al. 
(2002; hereafter B02) (dashed line) and Chavez et al (2012; 
hereafter Chl2) (solid line) are also shown. These data corre- 
spond to "young" single clusters in well studied giant extra- 
galactic Hll regions. They comprise the structureless entities 
in the regions free of bubbles and shells which are supposed 
to appear at more evolved phases of evolution. These lines 
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Figure 7. log(cr) — log(L) relations for the observed CNSFRs. Panels show the relations derived for the velocity dispersions of the gas 
using a single Gaussian fit (top), narrow components (centre) and broad components (bottom), respectively. The panels on the right 
show the result after correcting the luminosity for evolution (see text). The lines correspond to a linear fit to "young" Giant Hll regions, 
plotted as a reference value; B02 (dashed line) and Chavez et al. 2012 (solid line). 



may be taken as the locus of virialized structures sampled 
by these compact emission knots. 

Contrary to what is found for the compact emission 
line knots in nearby spirals by B02 and Chl2, no correlation 
is found between L(Ha) and a either for single Gaussian 
fits or for individual kinematical components; only three of 
our CNSFRs lie on top of the locus of virialized objects 
for the single Gaussian fit case. A weak correlation seems 
to emerge for the narrow kinematical component of hydro- 
gen that is maintained for the broad component although 
shifted to higher velocity dispersions. At any rate, the large 
dispersion of the data seems to be a common feature. 

A certain degree of scatter in the diagram is expected 
related to the age evolution of the ionizing regions which 



would lower their Ha luminosity while keeping their velocity 
dispersion almost constant. The presence of this effect is 
reinforced by the fact that the three regions on top of the 
B02 and Chl2 lines are those with the largest equivalent 
widths of H/3 (lower age) and therefore should be expected 
to be the less evolved ones. 

A correction for this effect can be made by using 
the same expression employed for the estimation of the 
mass of the ionizing clusters, using a reference value for 
the equivalent width of H/3 for an age lower than 3 Myr 
£;W"(H/3)o = 150A. 

r^TT„\ rCTJ„^ in[44.48 + 0.86 log[BW(ir/3)o] /~n 

L (Ha) Corrected = L(HajObsorvedlO' ""^ f-AU ^yj 

The panels on the right column of Fig 7 show the cor- 



Kinematical structure of CNSFRs 9 



relations just discussed once these evolutionary corrections 
have been applied. The net result is a considerable reduction 
of the scatter thus suggesting that the equivalent width of 
H/3 (i.e. age) is the second parameter underlying the L(Ha) 
vs. a relation. 

The stellar velocity dispersion behaves similarly to the 
broad component of the emission lines, albeit with larger 
dispersion. 

It is of major interest to find out how widespread is the 
presence of two distinct components in the emission lines of 
these ionised regions. Firstly, a change in position in the di- 
agnostic diagrams would certainly affect the classification of 
the activity in the central regions of the concerned galaxies. 
Secondly, it will affect the inferences about the nature of the 
source of ionisation in the two components. Thirdly, it could 
have an influence on the gas abundance determinations given 
that the ratio of the narrow to the broad components of H/3 
is, in most cases, about 1/2 . 

Clearly it is not possible to use global line ratios to 
estimate gaseous abundances if the permitted and forbidden 
line fluxes are originated in different kinematic systems. 

We analyse the relation between the Ha luminosi- 
ties of the broad component and the derived masses for 
the CNSFRs. Fig. 8 shows the relation between the log- 
arithm of the upper limit to the dynamical masses and 
log(I/(Ha)(,), together with the dynamical (Mdyn) and Kep- 
lerian (Mkbp) masses of Wolf-Rayet (WR) galaxies derived 
by (Lopez-Sanchez 2010, hereafter LSIO). Values of Mdyn 
were estimated from Hi radio observations considering the 
inclination-corrected maximum rotation velocity which is 
obtained at the maximum radius observed in their deep opti- 
cal images, and assuming virial equilibrium. As was pointed 
out by LSIO, since the extension of the neutral gas is usually 
larger than the extension of the stellar component, their val- 
ues of Mdyn may be underestimated, although they can be 
used as a rough estimation of the total mass of the galax- 
ies. The Mkbp were estimated from the kinematics of the 
ionised gas. The Keplerian masses are lower than the dy- 
namical ones, as expected, in almost all cases (see Fig. 11 
of LSIO), with the more massive galaxies showing higher 
MkhpI Mdyn ratios. This fact indicates that the kinematics 
derived from the ionised gas is not completely appropriate 
for deriving the total dynamical mass in this kind of objects 
(LSIO). 

The upper limits to the dynamical masses of the CNSF 
complexes seem to follow a sequence when we compare them 
with those masses derived for the WR galaxies (see Fig. 8). 
Our objects are located in the lower part of the sequence, 
showing lower masses and lower Hq luminosities. A fitting to 
the logarithmic dynamical mass-luminosity (M-L) relation 
shown by the CNSFRs (solid line in the figure) gives the 
following expression: 



log(L(Ha)6) = (0.8 ± 0.1) log(M.) + (33 ± 1). 



(8) 



The relation between the masses of the ionising clusters 
of the CNSFRs and their Hq luminosities is shown in Fig. 9. 
In spite of depending on a theoretical relationship which is 
only function of the Hq luminosity and the equivalent width 
of the H/? emission line we can appreciate that the CNS- 
FRs (including the objects studied by D07) seem to define 
a different sequence to the one followed by the Hll galaxies 
studied in Hagele et al. (2006; 2008; 2011, hereafter H06, 



H08 and Hll), and the WR galaxies analysed by LSIO (see 
Lopez-Sanchez & Esteban 2008, for a detailed description 
of the optical data) . This could be due to underestimating 
the EW(H/3) hence overestimating the masses of the ionis- 
ing star clusters. The dilution of the EW(H/3) can be caused 
by the presence of an old underlying stellar population with 
wider H/3 absorption that rises the continuum and depresses 
the emission line (see a detailed discussion about this effect 
in D07). In addition, the low EW(H/3) values for CNSFRs 
might come from the contribution to the observed contin- 
uum by the galactic bulge population. Dors et al. (2008) 
suggest that the combination of an underlying older popu- 
lation from the region itself and stars from the bulge could 
account for the observed effect. The strong influence that the 
measured equivalent width has on the derived mass of the 
ionising cluster indicates that the trend observed in Figs. 8, 
9 and 10 can be attributed to the fact that the contamina- 
tion by the bulge will be relatively stronger as the intrinsic 
luminosity of the starburst diminishes. 

This effect is more evident for the less massive objects. 
The most massive CNSFRs tend to be located around the 
position occupied by the Hll and WR galaxies of similar 
cluster ionising masses. The effect of the presence of an un- 
derlying old stellar population over the equivalent width of 
H/3 in this kind of emission line galaxies with recent episodes 
of star formations is relatively small since the continuum of 
the bright young stellar population dominates the light in 
this spectral range [see a detailed discussion about the con- 
tributions of the different stellar populations to the EW(H/3) 
in Perez-Montero et al. 2010 and Hll]. 

The solid line in Fig. 9 represents the fit to the loga- 
rithmic relation between Mton and total Hq luminosities for 
the CNSFRs, and is given by 

log(L(HQ)t) = (1.5 ± 0.4) log(M,„„) -f (30 ± 2). (9) 

Finally, in Fig. 10 we have plotted the relation between 
the masses of the ionised gas and the total Hq luminosities, 
together with the CNSFRs presented in D07 and the Hll 
galaxies studied in H06, H08 and Hll. The fitting to this 
relation for the CNSFRs gives 

log(L(HQ)t) = (0.67 ± 0.07) log(A/Hii) + (36.4 ± 0.3) (10) 



4 SUMMARY AND CONCLUSIONS 

We have analysed the results derived by H07, H09 and HIO 
using high spectral resolution data of 17 CNSFRs belonging 
to three barred spiral galaxies: NGC 2903, NGC3310 and 
NGC3351. The stellar velocity dispersions were measured 
from the CaT lines at AA8494, 8542, 8662 A. The gas veloc- 
ity dispersions have been measured by Gaussian fits to the 
H/3 A 4861 A and the [Qui] A5007A emission lines. 

Stellar velocity dispersions range between 31 and 
80kms~\ In the case of NGC 2903 and NGC 3351 (a. be- 
tween 37 and 64, and between 39 and 66kms~^, respec- 
tively) these values are about 25kms~^ larger than those 
measured for the gas from the H/3 emission line using a sin- 
gle Gaussian fit. For NGC 3310 the stellar (between 31 and 
80kms^^) and gas (H/?) velocity dispersions are in relatively 
good agreement, the stellar ones being marginally larger. 

Single Gaussian fitting to the [Om] 5007 A line provides 
gas velocity dispersions for NGC 2903 and NGC 3351 that 
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Figure 8. Ho luminosity vs. cluster dynamical mass for the observed CNSFRs (squares). Similar data on WR galaxies from the sample 
of LSIO are shown for comparison. In this latter case both dynamical (inverted triangles) and Keplerian (circles) masses are plotted. 
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are very similar to - or marginally larger than - the stel- 
lar ones. The gas velocity dispersion derived from a single 
Gaussian fit for the [Olll] line in NGC 3310 is also compara- 
ble to the stellar velocity dispersion but, contrary to what 
was found in the other two galaxies, its value is very similar 
also to that derived using H/?. 

The best Gaussian fits for H/3, however, involved two 
different components for the gas: a "broad component" with 
a velocity dispersion similar to that measured for the stars 
for NGC 2903 and NGC 3351, and about 20kms"^ larger 
than the stellar one for NGC 3310, and a "narrow compo- 
nent" with a velocity dispersion lower than the stellar one 
by about 30kms^^. This narrow component shows a rel- 
atively constant value for the two emission lines (H/3 and 
[Olll]), close to 23kms^^ (with scatter errors of 1.1 and 
1.4 kms"'^, respectively) for all the studied CNSFRs. 

The [Olll]/H/3 ratio distribution shows that the two 
systems are clearly segregated for the high-metallicity re- 
gions of NGC 2903 and NGC 3351, with the narrow compo- 
nent having the lowest excitation. In the regions of the low- 
metallicity galaxy, NGC 3310, these two components and 
those values derived using the single Gaussian fit are very 
similar. 

Values for the upper limits to the dynamical masses es- 
timated from the stellar velocity dispersion using the virial 
theorem for the studied CNSF complexes are in the range 
between 4.9 x 10® and 1.9 x 10* M©. Masses of the ionising 
stellar clusters of the CNSFRs have been derived from their 
Ha luminosities and the equivalent width of the H/3 emis- 
sion line under the assumption that the regions are ionisa- 
tion bound and have a single stellar population, and without 
taking into account any photon absorption by dust. These 



masses of the regions studied in the three galaxies vary be- 
tween 8.0 X lO'' and 4.9 x 10** M©. Therefore, the ratio of the 
ionising stellar population to the total dynamical mass, un- 
der these hypothesis, is between 0.01 and 0.16. The derived 
masses for the ionised gas, also from their Ha luminosities, 
vary between 7.0 x 10^ and 7.2 x 10^ Mq. 

The distribution of the CNSFRs in the log(o-) - log(L) 
plane presents a correlation between their luminosities and 
velocity dispersions, albeit with large dispersion. The single 
Gaussian fit for our CNSFRs are slightly shifted to lower lu- 
minosities and/or higher values of the gas velocity dispersion 
with respect to the regression found for virialised systems, 
although the general behaviour is similar to what was found 
by Firpo et al. (2010, 2011) for giant extragalactic Hll re- 
gions of NGC 6070 and NGC 7479, and star-forming regions 
belonging to the Haro 15 galaxy, respectively. The narrow 
components of the multi-Gaussian fits to the emission lines 
of the CNSFRs lie very close to this linear regression, with 
the brightest objects around the regression, while the broad 
components are located in a parallel sequence shifted to- 
wards lower luminosities and/or higher velocity dispersions. 
We have been able to considerably reduce the observed scat- 
ter in the L-cr relation applying an evolutionary correction 
that depends on the age of the regions through their equiva- 
lent widths of H/3. The EW(H/3) seems to be a second order 
parameter underlying this relation. Although the stellar ve- 
locity dispersions have a large scatter, they present a similar 
behaviour to the broad component ones. 

We have also analysed the relations between the total 
Ha luminosities and the derived masses for the studied CNS- 
FRs. The upper limits to the dynamical masses of the CNSF 
complexes seem to be located in the lower part of a sequence 
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showing lower masses and lower Ha luminosities than those 
masses derived for WR galaxies by LSIO. The distribution 
of the masses of the ionising clusters of the CNSFRs, includ- 
ing the regions presented by D07, falls outside a sequence 
defined by the Hll and WR galaxies analysed by H06, H08 
and Hll, and by LSIO, respectively. The CNSFRs are also 
arranged in a steeper slope sequence which is located toward 
lower values of the Ha luminosities or higher masses, except 
for the most massive regions. This occurs although the ionis- 
ing stellar masses depend on a theoretical relationship which 
is a function only of the Ha luminosity and the equivalent 
width of the H/? emission line. Underestimating EW(H/3) 
would produce an overestimation of these masses. We have 
to remark that, according to our findings, the super star- 
clusters in CNSFRs seem to contain composite stellar pop- 
ulations (see e.g. H07 and D07). The contribution of stars 
from the bulge projected along the line of sight (Dors et al. 
2008) has also to be considered. Although the youngest pop- 
ulation of star-forming complexes dominates the UV light 
and is responsible for the gas ionisation, it represents only 
about 10 per cent of the total mass belonging to the region. 
This can explain the low EWs of the emission lines gener- 
ally measured in these regions. This may well apply to the 
case of other super star-clusters and therefore conclusions 
drawn from fits of single stellar population models should 
be taken with caution (e.g. McCrady et al. 2003; Larsen 
et al. 2004). On the other hand, the relation between the 
derived masses of the ionised gas of the CNSFRs and the 
Ha luminosities follow the same sequence defined by the Hll 
galaxies presented by H06, HOB and Hll. 

A possible scenario for understanding the behaviour 
of CNSFRs in the L-a and (Jgas-o'* diagrams is to assume 
that the narrow component of hydrogen belongs to an in- 
ner rotating disc, therefore its velocity dispersion is approx- 
imately constant for all the CNSFRs. On the other hand, 
the broad component of hydrogen would be related to the 
CNSF complexes themselves, and therefore its velocity dis- 
persion would be expected to be similar to the stellar one. 
(The velocity dispersion of coUisionally excited [Olll] is in 
some cases larger than this, a behaviour also seen in Seyfert 
2 galaxies, but not in Seyfert 1 (see Jimenez-Benito et al. 
2000). In that case, all the Ha derived quantities for the 
CNSFRs should use the luminosity of the broad component 
which is about one half of the total. The fact that there is 
a distinctive behaviour in the Z/(Ha) vs. Mion in disk Hll 
regions and CNSFRs in the sense of showing a larger mass 
for a given Ha luminosity could be related to the Hll regions 
being matter bounded, that is photons are escaping the re- 
gion. Incidentally (or not so) these photons would ionise the 
hydrogen in the inner disc and therefore would in part be 
responsible for the narrow component of hydrogen. All this 
would be consistent with the picture of some CNSFRs un- 
dergoing "residual" star formation that involves about 10% 
of the total dynamical mass. This behaviour would be re- 
lated to the evolution of the region with the ones showing 
the larger EW(H/3), that is the younger ones, having more 
hydrogen available and thus being less prone to photon es- 
cape. The SFRs for the CNSFRs should also be derived from 
the luminosity of the broad component of Ha and would be 
about half the ones derived from the total Ha fiuxes. This 
would explain the lack of detection at radio wavelengths 
(Hagele et al. 2010). 



The existence of more than one velocity component in 
the ionised gas corresponding to kinematically distinct sys- 
tems, deserves further study. Several results derived from 
the observations of the different emission lines could be af- 
fected, among others: the classification of the activity in the 
central regions of galaxies, the inferences about the nature of 
the source of ionisation, the gas abundance determinations, 
the number of ionising photons from a given region and any 
quantity derived from them, as well as the cr-L and M-L re- 
lations. To disentangle the origin of these two components 
and to discriminate among the different stellar population 
contributions to the continuum of the regions it will be nec- 
essary to map these regions with high spectral and spatial 
resolution and much better S/N ratio. High resolution 2D 
spectroscopy with IFUs would be the ideal tool to approach 
this issue. 



ACKNOWLEDGEMENTS 

We are grateful to the referee. Angel Lopez-Sanchez, for a 
thorough reading of the manuscript and for his constructive 
comments and suggestions. 

The WHT is operated in the island of La Palma by the 
Isaac Newton Group in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofisica de Ca- 
narias. We thank the Spanish allocation committee (CAT) 
for awarding observing time. 

Some of the data presented in this paper were obtained 
from the Multimission Archive at the Space Telescope Sci- 
ence Institute (MAST). STScI is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555. Support for MAST for non- 
HST data is provided by the NASA Office of Space Science 
via grant NAG5-7584 and by other grants and contracts. 

This research has made use of the NASA/IPAC Ex- 
tragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration and of the SIMBAD database, operated at 
CDS, Strasbourg, France. 

Financial support for this work has been provided by 
the Spanish Ministerio de Educacion y Ciencia (AYA2007- 
67965-C03-03 and AYA2010-21887-C04-03). Partial support 
from the Comunidad de Madrid under grant S2009/ESP- 
1496 (ASTROMADRID) is acknowledged. ET and RT are 
grateful to the Mexican Research Council (CONACYT) for 
support under grants CB2005-01-49847F and CB2008-01- 
103365. 



REFERENCES 

Amori'n R., Vi'lchez J. M., Hagele G F., Firpo V., Perez- 

Montero E., Papaderos P., 2012, Astrophys. J. Letters, 

754, L22 
Boily C. M., Langon A., Deiters S., Heggie D. C, 2005, 

Astrophys. J. Letters, 620, L27 
Bosch G, Terlevich E., Terlevich R., 2002, Mon. Not. R. 

Astron. Soc, 329, 481 (B02) 
Bottinelli L., Gouguenheim L., Paturel G, de Vaucouleurs 

G, 1984, Astron. and Astrophys. Suppl. Sries, 56, 381 



Kinematical structure of CNSFRs 13 



Cervino M., Vails- Gabaud D., Luridiana V., Mas-Hesse 
J. M., 2002, Astron. Astrophys., 381, 51 

Chavez R., Terlevich E., Terlevich R., Plionis M., Bresolin 
P., Basilakos S., Melnick J., 2012, Mon. Not. R. Astron. 
Soc, 425, L56 

Colina L., Gonzalez-Delgado R., Mas-Hesse J. M., Lei- 
therer C, 2002, Astrophys. J., 579, 545 

Groom S. M., Lawrence J. S., Bland-Hawthorn J., Bryant 
J. J., Fogarty L., Richards S., Goodwin M., and 18 co- 
authors, 2012, Mon. Not. R. Astron. Soc, 421, 872 

Gumming R. J., Fathi K., Ostlin G., Marquart T., Marquez 
I., Masegosa J., Bergvall N., Amram P., 2008, Astron. 
Astrophys., 479, 725 

de Vaucouleurs G., de Vaucouleurs A., Gorwin Jr. H. G., 
Buta R. J., Paturel G., Fouque P., 1991, Third Reference 
Gatalogue of Bright Galaxies. Volume 1-3, XH, 2069 pp. 7 
figs.. Springer- Verlag Berlin Heidelberg New York 

Diaz A. I., Alvarez-Alvarez M., Terlevich E., Terlevich R., 
Portal M. S., Aretxaga I., 2000, Mon. Not. R. Astron. 
Soc, 311, 120 

Diaz A. I., Gastellanos M., Terlevich E., Luisa Garcia- 
Vargas M., 2000, Mon. Not. R. Astron. Soc, 318, 462 

Diaz A. I., Terlevich E., Gastellanos M., Hagele G. F., 2007, 
Mon. Not. R. Astron. Soc, 382, 251 

Dors Jr. O. L., Storchi-Bergmann T., Riffel R. A., Schimdt 
A. A., 2008, Astron. Astrophys., 482, 59 

Esteban G., Vilchez J. M., 1992, Astrophys. J., 390, 536 

Firpo v., Bosch G., Hagele G. F., Diaz A. I., Morrell N., 
2011, Mon. Not. R. Astron. Soc, 414, 3288 

Firpo v., Bosch G., Hagele G. F., Diaz R., 2013, Astrophys. 
J. (submitted) 

Firpo v., Bosch G., Hagele G. F., Morrell N., 2010, Mon. 
Not. R. Astron. Soc, 406, 1094 

Fleck J.-J., Boily G. M., Langon A., Deiters S., 2006, Mon. 
Not. R. Astron. Soc, 369, 1392 

Garcia-Benito R., Diaz A., Hagele G. F., Perez-Montero E., 
Lopez J., Vilchez J. M., Perez E., Terlevich E., Terlevich 
R., Rosa-Gonzalez D., 2010, Mon. Not. R. Astron. Soc, 
408, 2234 

Garcia- Vargas M. L., Bressan A., Diaz A. I., 1995, Astron. 
and Astrophys. Suppl. Sries, 112, 35 

Gonzalez-Delgado R. M., Heckman T., Leitherer G., 
Meurer G., Krolik J., Wilson A. S., Kinney A., Koratkar 
A., 1998, Astrophys. J., 505, 174 

Graham J. A., Phelps R. L., Freedman W. L., Saha A., 
Ferrarese L., Stetson P. B., Madore B. F., Silbermann 
N. A., and 15 coauthors, 1997, Astrophys. J., 477, 535 

Hagele G. F., 2008, PhD thesis, Universidad Autonoma de 
Madrid 

Hagele G. F., Ascasibar Y., Richards A. M. S., Gardaci 
M. v., Vasquez J., Diaz A. L, Rosa Gonzalez D., Terlevich 
R., Terlevich E., 2010, Mon. Not. R. Astron. Soc, 406, 
1675 

Hagele G. F., Diaz A. L, Gardaci M. V., Terlevich E., Ter- 
levich R., 2007, Mon. Not. R. Astron. Soc, 378, 163 (H07) 

Hagele G. F., Diaz A. L, Gardaci M. V., Terlevich E., Ter- 
levich R., 2009, Mon. Not. R. Astron. Soc, 396, 2295 
(H09) 

Hagele G. F., Diaz A. L, Gardaci M. V., Terlevich E., Ter- 
levich R., 2010, Mon. Not. R. Astron. Soc, 402, 1005 
(HIO) 

Hagele G. F., Diaz A. I., Terlevich E., Terlevich R., Perez- 



Montero E., Gardaci M. V., 2008, Mon. Not. R. Astron. 
Soc, 383, 209 (H08) 

Hagele G. F., Firpo V., Bosch G., Diaz A. L, Morrell N., 
2012, Mon. Not. R. Astron. Soc, 422, 3475 

Hagele G. F., Garcia-Benito R., Perez-Montero E., Diaz 
A. I., Gardaci M. V., Firpo V., Terlevich E., Terlevich R., 
2011, Mon. Not. R. Astron. Soc, 414, 272 

Hagele G. F., Perez-Montero E., Diaz A. I., Terlevich E., 
Terlevich R., 2006, Mon. Not. R. Astron. Soc, 372, 293 
(H06) 

Harper Jr. D. A., Low F. J., 1973, Astrophys. J. Letters, 
182, L89-1- 

Ho L. G., Filippenko A. V., 1996a, Astrophys. J. Letters, 
466, L83-h 

Ho L. G., Filippenko A. V., 1996b, Astrophys. J., 472, 600 

James B. L., Tsamis Y. G., Barlow M. J., Walsh J. R., 
Westmoquette M. S., 2013, Mon. Not. R. Astron. Soc, 
428, 86 

James B. L., Tsamis Y. G., Walsh J. R., Barlow M. J., 
Westmoquette M. S., 2013, Mon. Not. R. Astron. Soc, 
430, 2097 

Jimenez-Benito L., Diaz A. L, Terlevich R., Terlevich E., 
2000, Mon. Not. R. Astron. Soc, 317, 907 

Kehrig G., Monreal-Ibero A., Papaderos P., Vilchez J. M., 
Gomes J. M., Masegosa J., Sanchez S. F., Lehnert M. D., 
and 13 coauthors, 2012, Astron. Astrophys., 540, All 

Kennicutt Jr. R. G., Keel W. G., Blaha G. A., 1989, Astron. 
J., 97, 1022 

Larsen S. S., Brodie J. P., Hunter D. A., 2004, Astron. J., 
128, 2295 

Leitherer G., Heckman T. M., 1995, Astrophys. J., Suppl. 
Ser., 96, 9 

Lopez-Hernandez J., Terlevich E., Terlevich R., Rosa- 
Gonzalez D., Diaz A., Garcia-Benito R., Vilchez J., Hagele 
G., 2013, Mon. Not. R. Astron. Soc, 430, 472 

Lopez-Sanchez A. R., 2010, Astron. Astrophys., 521, A63H- 
(LSIO) 

Lopez-Sanchez A. R., Esteban G., 2008, Astron. Astro- 
phys., 491, 131 

Lopez-Sanchez A. R., Esteban G., 2009, Astron. Astro- 
phys., 508, 615 

Lopez-Sanchez A. R., Esteban G., 2010a, Astron. Astro- 
phys., 516, A104-h 

Lopez-Sanchez A. R., Esteban G., 2010b, Astron. Astro- 
phys., 517, A85-f 

Lopez-Sanchez A. R., Esteban G., Garci'a-Rojas J., Peim- 
bert M., Rodriguez M., 2007, Astrophys. J., 656, 168 

Lopez-Sanchez A. R., Mesa-Delgado A., Lopez-Martin L., 
Esteban G., 2011, Mon. Not. R. Astron. Soc, 411, 2076 

Macchetto F., Golina L., Golombek D., Ferryman M. A. G., 
di Serego Alighieri S., 1990, Astrophys. J., 356, 389 

McGrady N., Gilbert A. M., Graham J. R., 2003, Astro- 
phys. J., 596, 240 

Miller J. S., Mathews W. G., 1972, Astrophys. J., 172, 593 

Monreal-Ibero A., Relano M., Kehrig G., Perez-Montero E., 
Vilchez J. M., Kelz A., Roth M. M., Streicher O., 2011, 
Mon. Not. R. Astron. Soc, 413, 2242 

Monreal-Ibero A., Vilchez J. M., Walsh J. R., Muiioz- 
Tunon G., 2010, Astron. Astrophys., 517, A27-I- 

Morgan W. W., 1958, Publ. Astron. Soc. Pac, 70, 364 

Pastoriza M. G., Dottori H. A., Terlevich E., Terlevich R., 
Diaz A. I., 1993, Mon. Not. R. Astron. Soc, 260, 177 



14 G. F. Hagele et al. 

Perez-Montero E., Garcia-Benito R., Hagele G. P., Diaz 
A. I., 2010, Mon. Not. R. Astron. Soc, 404, 2037 

Perez-Montero E., Vilchez J. M., Cedres B., Hagele G. P., 
MoUa M., Kehrig C, Diaz A. I., Garcia-Benito R., Martin- 
Gordon D., 2011, Astron. Astrophys., 532, A141-I- 

Perez-Olea D., 1996, PhD Thesis. Universidad Autonoma 
de Madrid 

Pizzella A., Corsini E. M., Vega Beltran J. C., Bertola P., 
2004, Astron. Astrophys., 424, 447 

Planesas P., Colina L., Perez-Olea D., 1997, Astron. Astro- 
phys., 325, 81 

Relaiio M., Beckman J. E., 2005, Astron. Astrophys., 430, 
911 

Rieke G. H., Lebofsky M. J., 1978, Astrophys. J. Letters, 
220, L37 

Rieke G. H., Low P. J., 1972, Astrophys. J. Letters, 176, 
L95-h 

Rosales-Ortega P. P., Diaz A. I., Kennicutt R. C., Sanchez 
S. P., 2011, Mon. Not. R. Astron. Soc, 415, 2439 

Rosales-Ortega P. P., Kennicutt R. C, Sanchez S. P., Diaz 
A. I., PasquaU A., Johnson B. D., Hao C. N., 2010, Mon. 
Not. R. Astron. Soc, 405, 735 

Salpeter E. E., 1955, Astrophys. J., 121, 161 

Sanchez S. P., Kennicutt R. C., Gil de Paz A., van de Ven 
G., Vilchez J. M., Wisotzki L., Walcher C. J., Mast D., 
and 64 coauthors, 2012, Astron. Astrophys., 538, A8 

Sanchez S. P., Rosales-Ortega P. P., Kennicutt R. C, John- 
son B. D., Diaz A. L, Pasquah A., Hao C. N., 2011, Mon. 
Not. R. Astron. Soc, 410, 313 

Sanchez S. P., Rosales-Ortega P. P., Marino R. A., Iglesias- 
Paramo J., Vilchez J. M., Kennicutt R. C., Diaz A. L, and 
13 coauthors, 2012, Astron. Astrophys., 546, A2 

Sersic J. L., Pastoriza M., 1965, Publ. Astron. Soc. Pac, 
77, 287 

Sersic J. L., Pastoriza M., 1967, Publ. Astron. Soc. Pac, 
79, 152 

Telesco C. M., Harper D. A., 1980, Astrophys. J., 235, 392 



